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ABSTRACT 

Most of the observational studies of supernova (SN) explosions are limited to early phases (< a few yr 
after the explosion) of extragalactic SNe and observations of SN remnants (> 100 yr) in our Galaxy 
or very nearby galaxies. SNe at the epoch between these two, which we call "transitional" phase, 
have not been explored in detail except for several extragalactic SNe including SN 1987A in the Large 
Magellanic Cloud. We present theoretical predictions for the infrared (IR) dust emissions by several 
mechanisms; emission from dust formed in the SN ejecta, light echo by circumstellar and interstellar 
dust, and emission from shocked circumstellar dust. We search for IR emission from 6 core-collapse 
SNe at the transitional phase in the nearby galaxies NGC 1313, NGC 6946, and MlOl by using the 
data taken with the AKARI satellite and Spitzer. Among 6 targets, we detect the emission from SN 
1978K in NGC 1313. SN 1978K is associated with 1.3 x IQ-^Mq of silicate dust. We show that, among 
several mechanisms, the shocked circumstellar dust is the most probable emission source to explain the 
IR emission observed for CSM-rich SN 1978K. IR emission from the other 5 objects is not detected. 
Our current observations are sensitive to IR luminosity of > lO'^^erg s^^, and the non-detection of 
SN 1962M excludes the existence of the shocked circumstellar dust for a high gas mass-loss rate of 
~ 10^"* Mq yr~^. Observations of SNe at the transitional phase with future IR satellites will fill the 
gap of IR observations of SNe with the age of 10 — 100 years, and give a new opportunity to study the 
circumstellar and interstellar environments of the progenitor, and possibly dust formation in SNe. 
Keywords: dust, extinction — supernovae: general — supernovae: individual (SNe 1909A, 1917A, 
1951H, 1962M, 1968D, 1978K) — infrared: stars 



1. INTRODUCTION 

Core-collapse supernovae (SNe) are the explosions of 
massive stars at the end of their lives. SN ejecta ex- 
pand into the interstellar medium (ISM) or circumstel- 
lar medium (GSM) with a huge kinetic energy (~ 10^^ 
erg). First, SN ejecta experience a free expansion. As the 
ejecta sweep up the surrounding material, the SN begins 
to decelerate. When the swept-up mass becomes compa- 
rable to the mass of the SN ejecta at t ^ a, few hundreds 
yr (hereafter t denotes the time after the explosion) , the 
SN ejecta expand self-similarly (so-called Scdov phase). 
Finally, the SN merges into the ISM at t > 10, 000 yr. 

Observational studies of SNe have been performed 
mainly by the two distinct ways. One is the observa- 
tions of extragalactic SNe at their early phases {t < 10 
yr) , and the other is observations of supernova remnants 
(SNRs, t > 100 yr) in our Galaxy or very nearby galax- 
ies. In contrast, SNe at the epoch between these two 
(i = 10 — 100 yr), which we call "transitional" phase, 
have not been explored in detail. At such epochs, only 
a few objects have been studied. The best case is SN 
1987A in the Large Magellanic Cloud (|McCravl l2007l 
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and references the rein) . Other examples i nclude several 
X-ray SNe (e.g.. Ilmmler fc Kuntd 120051: llmmler et all 
[200a iSoria fc Pernal I2008D . and long-lasting SNe , such 
as SNe 1978K (IRvder etall IT993t ISchlegel et all p 99l 
[2000l [200llLenz fc Schleg elll2007':IS mith et al.ll2007ll and 



1988Z (jWiUiams et al.l l2002: Sc.hlegcl TTetrdl2006ll . See 



Sugerman et al.l 



also a recent paper on SN 1980K by 
(|20T1) . 

This transitional phase is important to understand a 
long ter m evolution of SNe, f rom SN to SNR. A recent 
study bv lLarsson et all ()2011l ) shows that SN 1987A, in 
fact, experiences a transition from SN to SNR in this 
phase. However, since most observations of SNe at tran- 
sitional phase except for SN 1987A have been performed 
at X-ray or radio wavelengths, the entire properties of 
SNe at the transitional phase are not fully understood. 

In this paper, we present the results of our search for 
infrared (IR) emission from core-collapse SNe at the tran- 
sitional phase. IR emission from SNe is expected to arise 
from dust associated with SNe. There are 3 possible 
populations of dust grains, according to their locations; 
(1) dust formed in the ejecta of SNe (SN dust), (2) dust 
formed by the pre-SN stellar wind, which is now located 
in the circumstellar region (CS dust), and (3) interstel- 
lar dust (IS dust). CS and IS dust can contribute to the 
IR emission by absorbing the radiation from SNe and 
reemitting it (light echo). CS dust is also heated and 
can also emit IR emission when it is swept up by the SN 
shock. In these ways, dust grains in and around SNe con- 
vert the large kinetic energy and luminosity of the SNe 
into the IR luminosity. Thus, IR observations of SNe can 
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Table 1 

List of Targets 



Host Galaxy Distance SN SN Type Position Position Rcf. 

(Mpc) a (J2000.0) & (J2000.0) 



NGC 1313 4.13 ±0.11 SN 1962M II 03hl8ml2.2s -66°31'38" 1 

SN 1978K Iln 03hl7m39.0s -66°33'048" 2 

NGC 6946 5.6 ±0.3 SN 1917A II 20h34m46.90s +60°07'29.08" 3 

SN 1968D II 20h34m58.41s +60°09'34.48" 4 

MlOl 6.7 ±0.3 SN 1909A'' II peculiar 14h02m03.1s +54°27'58" 3 

SN 1951H II 141i03m55.3s +54°21'41" 3 



References . — fl) lArtvukhina et al.|[T99g . (2 lDopita fc Hvdeil [T990I . (3) IBarbon et al.|[T999l . (4) 

Ivan Dvk et~aIlfT99l 

SN 1909A is out of field of view in the IRC images, and only in the Spitzer/MIPS image. 



Table 2 

Log of AKARI/IRC Observations 



Galaxy Date Observation ID Exposure Time (s) 

NIR, MIR-S MIR-L NIR, MIR-S MIR-L NIR MIR-S, L 



NGC 1313 2006 Dec 3 2006 Dec 4 1400416-001 1400417-001 178 196 
NGC 6946 2006 Dec 18 2007 Jun 16 1400620-001 1402217-001 178 196 
MlOl 2007 Jun 14 2007 Jun 17 1402211-001 1402212-001 178 196 



be used as a probe of the progenitor environment. 

In Section [2j we describe the target selection, the data 
used in the analysis, and results of photometry. We use 
the data taken with tw o IR satelli tes, AKARI satellit e 
([Murakami et al.ll2007| ) and Spitzer (| Werner et al.ll2004[ ). 
which provide the deepest images in near-IR (NIR) and 
mid-IR (MIR) wavelengths. Dust models used in this pa- 
per are shown in Section [3l Among 6 targets, we detect 
IR emission from SN 1978K. The origin of this emis- 
sion is discussed in Section H) In Section [5l theoretical 
predictions for the IR emission from normal SNe at the 
transitional phase are summarized, and implications of 
non-detection of 5 SNe are discussed. Prospects for fu- 
ture observations are discussed in Section |6l Finally, we 
summarize the conclusions in Section [T] 

2. TARGETS AND OBSERVATIONS 

2.1. Targets 

To search for IR emission from SNe at the transitional 
phase, we target 7 core-collapse SNe, namely SNe 1909A, 
1917A, 1951H, 1962M, 1968D, and 1978K, discovered in 
three very nearby galaxies NGC 1313, NGC 6946, and 
MlOl. These targets are selected based on the proximity 
to the host galaxy and the availa bility of NIR and M IR 
imaging data taken with th e IRC (Onaka ct al."2007D on- 
board the AKARI satellite (|Mura kami ct al. 2007). Gen- 
eral information on our targets is summarized in Table 

m 

Note that SN 1948B (Type II) in NGC 6946 is excluded 
because of the uncertainty in its position. SNe 1939C and 
1969P in NGC 6946 are also excluded since SN 1939C is 
classified as merely Type I (either core-collapse SNe or 
Type la SNe) and the type of SN 1969P is unknown. SN 
1980K is also located in NGC 6946. But since it is found 
th at AKARI data are not as deep as the data presented 
bv lSugerman et all (|2012D . we do not include SN 1980K. 
SN 1970G (Type IIL) in MlOl is excluded because the 
SN position is heavily contaminated by the H II region 



NGC 5455 (lFesenlll993D . 

For the distance to NGC 1313, we adopt 4.13 ± 0.11 
Mpc, which is estimated from the tip of the red giant 
branch jMcndcz ct al, 2002.') ■ We use t his value since it 
was al so used bv lLenz fc Schlegell (|2007| ) and lSmith et alJ 
(I2OOI for the study of SN 1978K. For NGC 6946, there 
are several distance measurements. We adopt the av- 
erage value of 5.6 zb 0.3 Mp c from the the H I Tully- 
Fishe r relation (lPiercelll994l) the CO TuUy-Fisher rela- 
tion (jSchoniger fc Sofuelll994| ). and the expandin g pho- 
tosph ere method of Type II SN 1980 K (jSchmidt et al.l 
IT99I and SN 2004et (|Sahu et al.ll2006[) . The distance to 
MlOl is assumed to be 6.7 ±0.3 Mpc, which is esti mated 
by the Cepheid calibration (jFreedman et al.ir2001[ ). 

2.2. AKARI IRC Data 

The AKARI/IRC equips three channels, i.e., NIR 
(1.8-5.5 Aim), MIR-S (4.6-13.4 /xm) and MIR-L 
(12.6—26.5 ^m). Each channel has about a 10' x 10' 
field of view. The NIR and MIR-S share the same field 
of view while the MIR-L observes a sky about 25' away. 
Thus, the data we use in this paper consist of two obser- 
vational sequences. The pixel scales of the detectors are 
1.46" X 1.46", 2.34" x 2.34", and 2.51" x 2.39" for the 
NIR, MIR-S and MIR-L, respectively 

We use the archived imaging data taken in the AKARI 
mission progr am " ISM in our Gala xy and Nearby galax- 
ies" (ISMGN: iKaneda et al]l2009( ). The data are taken 
with the two-filter mode [ Astronomical Observation 
Template (AOT) IRC02, see lOnaka et"all (12001 ]. We 
use 6 band images in total, i.e., N3 (reference wavelength 
3.2 ^m), N4 (4.1 /im), S7 (7.0 ^im), Sll (11.0 //m), L15 
(15.0 /im), and L24 (24.0 //m). The wavelength coverage 
with the response larger than 1/e of the peak is 2.7—3.8 
^m, 3.6-5.3 ^m, 5.9-8.4 ^m, 8.5-13.1 /xm, 12.6-19.4 
/xm, and 20.3-26.5 ^im for N3, N4, S7, Sll, L15, and 
L24, resp ectively. For more details of the IRC instru- 
ment, see lOnaka et al.l ()2007D . The log of observations is 
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Table 3 

Summary of Spitzcr/MIPS Data 



Galaxy 


Date 


Exposure Time (s) 


Observing mode 


AOR key 


Program ID 


NGC 1313 


2007 Sep 29 


42 


Medium Scan 


22618112 


40204: Kennicut et al. 


NGC 6946 


2007 Jul 10 


493 


Photometry 


18271232 


30494: Sugerman et al. 


MlOl 


2008 Jan 5 


42 


Medium Scan 


21380144 


40352: Rickc et al. 



Table 4 

Summary of Photometry 





SN 


Age ^ 






Flux (mjy)'^ 












(yr) 


IRC/N3 


IRC/N4 


IRC/S7 


IRC/Sll 


IRC/L15 


IRC/L24 


MIPS/24/im 


SN 


1962M 


44.0 
44.9 


<0.2 


<0.2 


<0.2 


<0.2 


<0.9 


<2 


<1 


SN 


1978K 


28.4 
29.2 


0.074 (0.014) 


0.086 (0.010) 


0.287 (0.039) 


2.85 (0.24) 


2.53 (0.13) 


2.73 (0.20) 


3.06 (0.40) 


SN 


1917A 


89.5 
90.0 
90.0 


<0.1 


<0.2 


<8 


<6 


<10 


<10 


<10 


SN 


1968D 


38.8 
39.3 
39.4 


<0.3 


<0.3 


<4 


<5 


<6 


<8 


<5 


SN 


1909A 


99.0 














<0.7 


SN 


1951H 


55.9 
56.3 


<0.07 


<0.1 


<0.4 


<0.6 


<1 


<3 


<i 



Age is given as year after the discovery, which is a sound approximation for the age after the explosion. 
^ Values in parenthesis for SN 1978K represents la error including the measurement error and the calibration error (5 %). 
In the case of non-detection, 3(t upper limit is given. 

summarized in Table [2] 

The data were reduced by using the IRC imaging 
pipel ine software version 20091022 (see iLorente et al.l 
120071 for details). The pipeline performs the basic re- 
duction, e.g., correction of bad pixels, subtraction of dark 
current, rejection of cosmic rays, linearity correction, flat 
fielding, and co-adding of the frames. 

To determine the accurate position, we match the 
sources in the N3, N4 and S7 band im ages with those 
in Tw o Micron All Sky Survey catalog (jSkrutskie et al.l 
|2006[) . The uncertainty in the position is < 2" for the N3 
and N4 band images, and < 3" for the S7 band images. 
Then, the Sll, L15 and L24 images are matched with S7 
images. In these images, the uncertainty in position is 
estimated as ^ 3", which is larger than that in the N3, 
N4 and S7 images because the matching is based on a 
small number of common sources in the field of view. 

To obtain the flux densities of the sources, we perform 
aperture photometry with the aperture of 10 pixel and 
7.5 pixel for the NIR and MIR-S/L bands, respectively. 
As the background, we use the annulus of 5 pixel width 
just outside of the aperture radius. When the SN is lo- 
cated at the edge of the image, the apertur e correction is 
applie d. The flux is calibrated by foUowing lLorente et al.l 
(j2007l) . Color correction is not performed to the derived 
photometry because the intrinsic spectral shape of the 
source is not known. Instead, we apply correction to the 
model when it is compared with the observations (see 
Section [3]) . 

When the SN is not detected, we derive an upper limit 
of the flux by putting an artificial source and measuring 
the flux of the source. For this purpose, we need a point 
spread function (PSF) of each image. For the N3 and N4 



images, PSF is constructed by the point sources in the 
same image. Since there are not many sources in the S7, 
Sll, L15 and L24 images, we use an av erage PSF con- 
struc ted from sources in different fields (jArimatsu et al.l 

EoliD. 

2.3. Spitzer MIPS Data 

In addition to the AKARI/IRC data, we use archive 
data taken with t he Multiband Ima ging Photometer 
for Spitzer (MIPS iRieke et al.l [200l onboard Spitzer 
(jWerner et al.l[200l . 

All the 6 SNe are within the field of view of some ob- 
servations with the MIPS 24 fim band. NGC 1313 was 
observed with the scan mode in the Spitzer Local Volume 
Legacy program (jPale et al.l '20091. MlOl was also ob- 
served with the scan mode bv Gordon et al.l (pOOg ). SNe 
1917A and 1968D in NGC 69 46 are within the field of 
SN 2002hh (jMeikle et al.ll2006f ). A summary of the data, 
including the Astronomical Observation Request (AOR) 
keys. Program IDs and PI names, is shown in Table [3l 

We use the post-basic calibrated data (PBCD) prod- 
ucts, which are produced by the Spitzer pipeline. The 
pixel scale of the final image is 2.45" x 2.45". The ref- 
erence wavelength of the MIPS 24 ^m band is 23.68 fim 
and the wavelength coverage is 20.5 fim - 28.5 /xm (the 
range where the response is > 10% of the peak). The 
positional uncertainty is 1.4". 

For photometry, we perform aperture photometry. 
When the SN is not detected, we derive an upper limit 
of the flux by putting an artificial source around the SN 
position. We constructed a PSF from point sources in 
the scan mode images of NGC 1313 and MlOl and use 
it to derive the upper limits. 



4 



Tanaka et al. 



< 



-30 -60 -90 




o 




o 








o 




o 














o 


o 

V 


o 


ro 


0) 
















o 


o 


o 


o 


u 


o 


"? 


o 


m 


o 


< 


1 

o 


CO 




(O 

1 


o 




o 


O) 




Ol 



1 90 60 30 -30-60-9C) 
A RA (arcsec) 



-30 -60 -90 
o 




1 90 60 30 -30 -60 -9C) 
A RA (arcsec) 



< 



-60 -90 
o 




1 90 60 



-30 -60 -9C) 



A RA (arcsec) 



30 -30 -60 -90 
o 




1 90 60 30 -30 -60 -9C) 
A RA (arcsec) 



< 



-30 -60 -90 
o 




90 



60 30 -30 -60 -9CI 
A RA (arcsec) 



o 
ro 

I 

o 

(O 

I ■ 
o 
a> 



1 90 



60 



30 



-30 

^«i-I-r 



-60 -90 
o 

O) 



J I I I I I I l-J I I I I I I I L. 



60 30 
A RA 



o 
I 

o 
<o 
I 

■ o 

O) 



-30 -60 -9C) 
(arcsec) 



Figure 1. AKARI/IRC images around SN 1978K in NGC 1313 (180" x 180" section, which corresponds to about 3.6 kpc X 3.6 kpc). 
North is up, and cast is to the left. Note that the PSF of the N3 a nd N4 bands are sli ghtly elongated in images on orbit due to a 
misalignment of the telescope and color aberrations in the NIR channel ijKaneda et al.|[2007l ^. 



The Spitzer/MIPS data are found to be slightly deeper 
than the AKARI/IRC data at 24 (ScctionEJ]). Thus, 
we use MIPS data for the flux of 24 /im. 

2.4. Results of Photometry 

We detect IR emission associated with SN 1978K in all 
the bands (Tabic E]). Figure [1] shows the AKARI/IRC 
images around SN 1978K (180" x 180" section). The 
source has an extremely red color with a flux ratio of 
Fy(S7)/i^^(Sll) = 0.1, which is out of a typical range 
in the sp iral arm and the inter-arm of the galaxy (~ 
0.8 - 1.4. lSakon et al.ll2007l ). The origin of this emission 
is discussed in Section |4l 

For the other 5 objects, wc do not detect a significant 
signal from SNe. Upper limits of the flux are shown in 
Table m 

3. DUST MODELS 

In the following sections, we compare the observations 
with various model spectra. We consider two types of 
carbon dust grains, i.e., amorphous carbon and graphite, 
and two kinds of silicate dust grains, i.e., "astronomi- 
cal silicate" and forsterite (Mg2Si04). We assume that 
grains are spherical. 

Thermal emission flux from the dust toward a SN can 
be expressed as follows; 



-Fa = M, 



dust 



47r . 



dust ) 



f 1 



V 47rd2 



(1) 

where Mjust is the total dust mass, p and a arc, respec- 
tively, the density and the radius of the dust grain, Qf^^ 
is the absorption efficiency factor, T^ust is the temper- 



ature of the dust, and d is the distance to the SN. We 
assume a single temperature for simplicity. In the equa- 
tion above, all of the IR emission from dust grains is 
assumed to be detected without being absorbed. 

We consider a single size of dust a = 0.1 fim. In reality, 
there is a distribution in the grain size, and if dust grains 
are swept up by the shock waves of the SN, their size 
distribution can change with time (jNozawa et al. 2007.). 
However, Qf'^/a is independent of the grain size if A ^ 
27ra, and thus, the assumption of the single-size dust does 
not affect the estimate of the dust mass. 

For the grain densities, we adopt p = 2.26, 2.26, 3.8, 
and 3.2 g cm^"^ for amorphous carbon, graphite, astro- 
nomical silicate, and forsterite, re spectively. T he absorp- 
tion efficiency factor i s taken from lEdohl (|1983f ) for amor- 
phous carbon. [Praind (l2003bD fo r graph ite and astronom- 
ical silicate, and lSemenov et all POOS') for forsterite. 

Since the observed flux in Table His not color-corrected 
due to the unknown intrinsic spectrum of the sources, 
we apply the color corrections to the model. Since the 
flux density of the AKARI/IRC data refers to the spec- 
trum F\ oc 1/A, the following correction is applied to the 
model spectrum F^°'^''\X); 



/ °'^°'(A)i?(A)dA 
J{K/X)R{X)dX ■ 



(2) 



Here, Xi is the reference wavelength of the band, and 
i?(A) is the response function (electron energy"^). We 
use the response functi o n measured in the l aboratory 
(see lOnaka et all 120071 : ILorente et all I2007D . The cal- 

http:/ /www. ir.isas.jaxa.jp/AKARI/Observation/ 
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Figure 2. The IR-radio SED of SN 1978K (filled circles) com- 
pared with that of Cas A (open squares) at the same distance (4.13 
Mpc). The red dashed line shows the power-law fit to the IR and 
radio data of SN 1978K while the blue dashed line shows the fit 
only to the radio data. The radio data of SN 1978K ar e taken from 



f th ct al. (2007). T he dat a of Cas A are taken from 
77), MezKcr et al . 1198^). ILiszt 



ro^) ! and lBarlow et al.l 1 12010^ 



Lucas! ilTgggI ') 



Baars et al 



Hines et al 



ibration of the Spitzer/MIPS data uses the blackbody 
spectrum with Tq = 10000 K as the reference. Thus, we 
apply the following correction; 



F- 



mips 



A 



exp(/ic/AiA:ro) ^ 1 
exp{hc/ XkTo) — 1 



,(3) 



RiX)dX 



We use the response function provided by the Spitzer 
Science Center Q. 

4. IR EMISSION FROM SN 1978K 

4.1. Spectral Fit 

Among our 7 targets, we detect IR emission only 
from SN 1978K. SN 1978K is an extraordinarily strong 
Type Iln SN, which is visible both in r adio and X- 
ray even at ~ 3 yr after the exp losion ('Ryde r et al 
1993: Petr e et alST9 M: Clm gai et alTl 995: Schlcge l et al 
1996: Montes ct al. 1997; Schlcgel ct al. 1999; Chu ct al 
1999: Schlcgel etjil. 1000, 2004; Lenz k Schlcgel 2007, 
Smith et al.l I2007D . Such long-lasting radio and X-ray 



emissions are clear evidence for the presence of the dense 
circumstella r medium, with which the SN ejecta has been 
interacting (IChevalieiill982l: IChevalier &: Franssonlll994l; 
IChugai et al.lll995D . 

Figure [5] shows the spectral energy distribution (SED) 
of SN 1978K from IR to radio wavelengths (filled circles). 
The SED of Cassiopeia A (Cas A) at the distance of NGC 
1313 is also shown for comparison (open squares). These 
two objects show similar overall SEDs from IR to radio 
wavelengths. 

Since SN 1978K is powerful in radio wavelengths 
(jSmith et al.l 120071) . we must take into account the in- 
fluence of the synchrotron em i ssion on the IR bands 
(see e.g., iMezger et all 119861; iRho et al.l l2003t iDwekl 

http:/ /ssc. spitzer. caltech.edu/mips/calibrationfilcs/ 
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Figure 3. The IR SED of SN 1978K (open black circle) com- 
pared with the model spectra. The red line shows the spectrum 
of the amorphous carbon with T(juat = 180 K and M^ust = 
6.8 X 10~^Mq while the blue line shows the spectrum of astro- 
nomical silicate with T(jugt = 230 K and Mjust = 1-3 x IO'^Mq. 
The crosses show the color-corrected model flux for each band. 
For both models, the synchrotron component {Fi/ oc i/-0 &0^ jg 
added. In the case of astronomical silicate, the color corrections 
are F^'^'^''^'^^(Ai)/Fm°dol(;)^.) ^ i q g^^ ^^^^ n ^j, ^he 
S7, Sll, L15, and MIPS 24 fim band, respectively. In the case of 
amorphous carbon, these are 1.1, 0.93, 0.95, and 1.0, respectively. 



I200I iHines et al.ll200l iBarlow et alllMol for the case 
of Cas A). The radio emission of SN 197 8K is weh fit- 
ted by a p ower-law spectrum oc (|Montes et al.l 
Hill [Smith ct al. 2007). When we fit the f lux of 4 ra- 
dio bands of SN 1978K (iSmith et al.ll2007D . we derive 
a = -0.67±0.05 (blue dashed line in FigureE]). The flux 
densities at the N3 and N4 bands are almost in line with 
the power-law spectrum. If we include the flux of these 
2 bands for the fitting, we obtain a = —0.60 ± 0.02 (red 
dashed line). Since these two values are consistent within 
the error, we assume that the emissions at the N3 and 
N4 bands are of synchrotron origin and use a = —0.60 
in the following discussions. 

The flux density of the S7, Sll, L15, and L24 (MIPS 24 
fim) bands is clearly above the synchrotron component. 
Therefore, we conclude that this excess is caused by the 
dust emission associated with SN 1978K. 

We fit the observed flux with the model spectra of the 
4 dust species with the dust temperature and mass as 
parameters. The quality of the fit is judged by for the 
4 observational points, i.e., the IRC S7, Sll, L15, and 
MIPS 24 /im bands. For each dust species, we choose 
the best model that gives the smallest x^- Since the de- 
pendence on the temperature and mass is simple, the fit 
is unique, i.e., there is no other solution with a different 
set of temperature and mass. 

Among the 4 dust species, the smallest is obtained 
for astronomical silicate with Tdust = 230 K and Mdust = 
1.3x IO-^Mq (blue line in Figure^]). The total IR fiux by 
the dust emission is Fir = 7.0 x lO^^'^ erg cm~^ s~^ and 
the total IR luminosity is Lm = 1.5 x lO^^erg s~^. Sili- 
cate dust grains are favored over carbon grains because 
the observation shows excess at the Sll band, which can 
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be attribu ted to the Si-0 stretching mode of silicate dust 
at 10 /xm (jPraine fc LeellT98llDrainell2003al) . As shown 
in the red hue in Figure[31 amorphous carbon dust grains 
are not able to explain the excess. Forsterite fits the ob- 
servations quite well, but gives a band feature at 10 /.tm 
slightly narrower than astronomical silicate. Hereafter, 
we simply call the dust species found in SN 1978K "sili- 
cate dust" . 

4.2. Origin of the IR emission 

SN 1978K is associated with 1.3 x IQ-^Mq of sili- 
cate dust. The derived temperature is relatively high 
(230 K). In this subsection, we discuss possible origins 
of the IR emission. Since SN 1978K is a CSM-rich 
SN, we first consider CS dust. C S dus t can radiate 
IR emission by light echo (Section I4.2.ip . Then, the 
echoing CS dust is located at the distance larger than 
D ^ ct/2 ^ lA X 10^3 (i/30 yr) cm, where the SN for- 
ward shock cannot reach. On the other hand, the CS 
dust swept up by the SN forward shock can be heated 
via the collisions with energetic electrons and emits ther- 
mal radiation (Section l4.2.2p . Then, the typical location 
is at - 5 X 10"(t/30 yr)(usN/5000 k m s~^) cm. An- 
other possibility is the SN dust (Section l4.2.3p . We show 
that, among these possibilities, the shocked CS dust is 
the most probable origin for the IR emission from SN 
1978K. 

4.2.1. Light Echo 

The first possibility for the IR emission is echo emis- 
sion of CS dust. The dust can absorb the ultra-violet 
or optical photons o f a SN, and radiate IR emission 
eode fc EvansI IT980I: Imlghtl IT980I: iDweli [1981 IT981 
IGraham fc Meikld 119861) . In fact, the presence of CS 
dust ech o has been suggested by early phase observations 
of SNe (IMeikle et all I2OO6I [2007t [Andrews et all IMll) 
and S N 1980K at the transitional phase |Sugerman et al.l 
[20TI . 

We calcula t e the luminosity of the light echo by fol- 
lowing ([Ull) . The gas mass-loss rate of the pro- 
genitor of SN 1978K is estimated to be about M 
10" ' 



tions (IRvder et al.lll993l: IChugai et al]|1995l: [Smith et al.l 
[2007[) . The duration of mass loss is assumed to be 
iwind = 10^ yr. This is a rather extreme assumption 
since the progenitor loses mass as large as 10 Mq prior 
to the explosion. We assume that the wind velocity is 
Wwind = 10 km s~^, which gives outer radius of CSM 
about i?2 = 3 X 10^* cm. The mass-loss rate is assumed 
to be constant, and thus, the CSM density has a power- 
law radial profile oc r~^. 

For the SN luminosity, we use an exponential form 
L = Lq 6xp(t/isN). We assume a normal peak luminos- 
ity Ln = 6 X 10^^ erg s~ ^ and isN = 25 days (jSahu et al.l 
|2006l : lKotak et al.ll2009[ ). The other parameters are set as 
follows; the gas-to-dust mass ratio / = 200, the grain ra- 
dius a = 0.1 /im, the mean absorption efficiency Q = 1, 
and the evaporation temperature of dust T = 1500 K. 
The SN outburst luminosity will make a dust cavity with 
a radius of i?i = 7 x 10^^ cm. 

The left panel of Figure 0] shows the computed IR lu- 
minosity of the CS dust echo (solid black line). Optical 
data of SN 1978K are also plotted for comparison under 
a simple assumption of zero bolometric correction (blue 
and red points). The IR echo luminosity at t = 10 — 100 
yr depends largely on the outer radius {R2) of the CSM 
since the luminosity fr om the CS e cho declines rapidly 
at t = 2R2/C ^ 1 yr pweklll983h . To illustrate this 
dependence, we also show the model with a longer mass- 
loss duration (dashed line for <wind = 10^ yr). With 
twind ^ 5 X 10^ yr, with which the progenitor loses more 
than 50 Mq, the IR echo does not drop dramatically at 
t <30 yr. But, even with this extreme case, the IR lumi- 
nosity is only an order of 10'^^ erg s~^ at t = 30 yr, which 
is much less than the observed IR luminosity. Thus, we 
conclude that the CS dust echo of SN outburst is unlikely 
to be the origin of the IR emission. 

We also calculate the IS dust echo by the same method 
as that for the CS dust echo. In fact, IR echo by 
the I S dust has been detected for several extragalactic 
SNe (IMeikle et a l. 2007; Kotak e t al. 200 9; Meiklc et al] 



2011 1 ) and Cas A (jKrause et al.l 120051 : iDwek fc Arendt 



2008( 1. The difference from the CS dust is that the den- 



Mq yr-1 from optical. X-ray and radio observa- ^i^^ to be uniform and that the much larger outer 
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Figure 5. The mass of shocked CS dust that survives at t = 30 
yr after the explosion as a function of mass-loss rate. The dashed, 
solid and dotted lines shows the model with ao = 0.3, 0.1, and 0.01 
/im, respectively. The red dashed line shows the model without 
dust destruction with ao = 0.1 /im. The other parameters are set 
as follows; vsn = 4000 km s~^, I'wind = 10 km s~^, and / = 200. 

radius is allowed. Here, the distribution of IS dust is ap- 
proximated to be a sphere with the radius of i?2 = 100 
pc (t = 2R2/C ^ 700 yr). The black solid and dashed 
lines in the right panel of Figure |4] show the luminosity 
evolution of the IS dust echo for the ISM gas density of 
n = 0.1 and 1 cm^'^, respectively. 

With the normal SN luminosity, the expected IR lu- 
minosity at i = 30 yr is not high enough to explain 
the observed luminosity of SN 1978K. To explain the 
observation, the SN had to be very luminous, Lq ~ 
6 X 10'*^ erg s~^, or the absolute optical magnitude of 
about —20 mag. It would imply that the brightness of 
SN 1978K was of the 8th magnitude. T his is in contrast 
to the suggestion bv lRvder et all (| 19931) that SN 1978K 
was subluminous. Although a very large luminosity of 
SN 1978K is not completely excluded given the sparse 
observations in 1978, we conclude that IS dust echo is 
also unlikely. 

4.2.2. Shocked CS dust 

The second source for the IR emission is the CS dust 
heated by the forward shock of the SN. The mid-IR emis- 
sion from the ring in S N 1987A is thoug ht to arise from 
this mechanism (e.g., iDwek et al.l 120081) . It is interest- 
ing to note that the emission properties of SN 1987A ring 
and SN 1978K look similar. The dust temperature in the 
SN 1987A ring is 165 - 185 K at ^ = 17.9 - 19.7 yr 



(jBouchet et al.|[200l iDwek et al.|[200l iSeok et al.|[2008[ ) , 

while that in SN 1978K is 230 K at t = 29 yr. In addition, 
the IR to X-ray flux ratio is of the order of unity in both 
objects (For SN 1 978K, Fx ^ 1 x lO'^^ erg cm'^ g-^, 
ISmith et al.l 120071 ). The similarity seems to support the 
shock-heated CS dust as an origin of IR emission from 
SN 1978K. 

We estimate the mass of the emitting CS dust swept by 
the forward shock. CS dust is subject to the destruction 
by the SN forward shock. The timescale of sputtering is 
given by 



-1 



^ 10« yr ( (-^ 
\ 1/im /VI cm 



, -1 



(4) 



Arendt|[T99^ . where Ug is the number density 
of shocked gas. Then, the evolution of the dust radius 
is estimated as a(t) = ao(l — t/r^^^^), where rgp^t = 
10^ yr(ao/l/.tm)(ng/l cm~'^)~^. Hereafter, the subscript 
denotes the value before the dust destruction. The gas 
mass-loss rate and the CS dust density is related as fol- 
lows; M = 47ri?2pwindt'wind = 47ri?2/pdust,oWwind , where 
Pwind is the mass density of the wind gas. We also define 
the mass density of dust before the destruction; pdust.o = 
(4/3)7rpaondust,o, where 7idust,o is the number density of 
dust before the destruction. The CS dust density after 
the destruction is given by pdust(i) = (4/3)7rndust,opa^ = 

(4/3)7rndust,opao(l - V'^sput)^ = Pdust,o(l - V^sput)^- 

Using these relations, the mass of the emitting, shocked 
CS dust is described as follows; 



Mcsdust(t) = / 'inR^pdust{t')vsNdt' 



47ri?2pdust,0l'SN(l - t' /T^p^tfdt' 



8 X IQ-^Mq 



wind 



'^wind 



10 km s 



10-5 Mq yr-i y V5000 km s-i 
f A"7/o(l-i'M'put)=^di' 



200 



30 yr 



(5) 



At the limit of no dust destruction, the integral in the last 
parentheses becomes t, and the shock-heated CS dust is 
proportional to the mass-loss rate. But r^p^^ is a function 
of the shocked gas density, and thus, the mass-loss rate. 
With a larger mass-loss rate, the sputtering timescale 
becomes shorter, and the integral in the last parentheses 
becomes smaller than t. 

Figure [5] shows the shocked (but surviving) mass of 
the CS dust at t = 30 yr as a function of mass-loss rate. 
The dashed, solid and dotted lines in black show the 
models with uq — 0.3, 0.1, and 0.01 /im, respectively. The 
red solid line is the model without dust destruction with 
ao = 0.1 ^m. Other parameters are assumed as follows; 
wgN = 4000 km s-\ Wwind = 10 km s~\ and / = 200. 
The SN shock velocity of SN 1978K is estimated to be 
vsN < 400 km s~^by the VL BI image, which marginally 
resolve it (|Smith et al.ll2007[ ). Thus we assume wsn = 
4000 km s-^. 

For small mass-loss rates, the destruction is not so ef- 
fective that the shocked CS dust mass is roughly propor- 
tional to the mass-loss rate. For large mass-loss rates, 
the surviving dust mass is not sensitive to the mass-loss 
rate. Because of the high gas density, the destruction 
is effective enough that the mass of dust swept up by 
the shock compensates for the mass destroyed by the 
shock 0. In the calculation above, we assume there is 
no dust-free cavity created by the SN outburst. Since 
most of surviving dust is located near the SN shock (at 
i? - 5 X 10i^(t/30 yr)(i>sN/5000 km s'^) cm), which 
is larger than the cavity size ('^ 7 x IQ-'^^ cm). Thus, 

* IFox et al.l | |20TO| . [20Tll) also derived the expression for the 
shock-heated dust mass, which is independent on the mass-loss 
rate. Since they assume a destruction-free shell, the mass of CS 
dust is larger than our calculation by a factor of 4. 
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the results presented here are not largely affected by the 
presence of the dust-free cavity. 

From the IR SED of SN 1978K, we estimate the mass 
of emitting dust to be 1.3 x 10~^Mq. The model of 
the shocked CS dust with M = lO"" Mq yr'^ and 
fSN = 4000 km s~^ is consistent with the observed dust 
mass if the initial dust radius is on = 0.3 /xrn . A sim ilar 
conclusion has been reached by iDwek et al.l ()2008() for 
SN 1987A. 

We also check the expected temperature for the 
shocked CS dust. With M = lO"'' Mq yr-\ the shocked 
gas density is about 6 x 10'^ cm^"^. For the electron tem- 
perature about 10^ K, the expected dust temperature is 
Tdust ^ 270 and 2 20 K for the dust r adius a = 0.1 and 0.3 
/im, respectively (jPwek et al.ll2008l ). Thus, the expected 
temperature of the shocked CS dust is also in good agree- 
ment with the observations. The surviving, shocked CS 
dust is distributed in a thin shell, whose inner radius is 
about 95 % of the radius at the forward shock. Using the 
mass absorption coefficient of dust for astronomical 
silicate at 24 /im (|Drainc 2003bh . the optical depth of 
dust is about 4 x 10~^(K,y/562 cm^ g~^), which is consis- 
tent with our assumption of optically thin dust (Section 
[3]). In summary, the IR emission from the shocked CS 
dust can naturally explain the observed IR emission from 
SN 1978K. 

4.2.3. SN dust 

The last possibility for the origin of the IR emis- 
sion is dust formed in SN ejecta. Theoretically, 
SNe are expected to form dust masses as large as 
Mdust = 0.1-1.0 Mq a t the inner layer of the SN 
eiecta (iKozasa et al.| [l989l. fl99lt ITo dini fc Ferrara|[200ll: 



Nozawa et all 120031 12008 



, 2010: Bianchi fc Schneided 

2007t ICherchneff fc Dwd3 12010[) . In fact, such a large 
amount of dust is required to explain the observed dust 
mass in high-redshift quasars, if the dust is formed only 
by SN e (iDwek et al.ll2007t IMeikle et al.ll2007l: IGah et al.l 

MM- 

Observationally the mass of dust formed by SNe 
is quite controversial. IR observations of young SNe 
(t < 1 — 2 yr) in external galaxies usually find only 
10^5 - IQ-^ Mq of dust (e.g., [W ooden et"all [l99l 



Ercolano et al.ll2007l: lElmhamdi et ah 2003: Meikle et al 
2007t iKotak et al.ll2009l: iSzalai et al.feoiOt IMeikle et al 



2011[) . On the other hand, observations of supernova 
remnants (SNRs) in our Galaxy and Large/Small Mag- 
ellanic Clouds detect a larger amount of dust. For exam- 
ple, by intensive observations of Cas A at IR to sub-mm 
wavelengths, it is estimated tha t ~ O.lMm of dust is 
formed in total (IRho e t al.l 120081: ISibthorpe et"aL| [20l 
IBarlow et all 120101) . Recentlv. iMatsuura et al.l (|2011l 
detected a large amount of cold dust in SN 1987A at 
t = 23.3 yr in far-IR and sub-mm wavelengths. The es- 
timated dust mass is as large as ^ 0.4 — O.TMq, and the 
dust temperature is ~ 20 K. 

To estimate the IR emission from SN dust, we first con- 
sider the total luminosity available, which gives an upper 
limit of the IR luminosity. At the transitional phase, one 
of the major heating sources of SN dust formed at the 
inner layers of the ejecta is the radio active decay of ^^Ti, 
which gives L ~ 10'^^ — lO^^erg (jFransson fc Kozmal 
120021) . This is much smaller than the IR luminosity 
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Figure 6. Expected IR luminosity for normal SNe. The red and 
blue lines show the echo of the IS dust and the CS dust, respec- 
tively. For the IS dust echo, the red solid and dashed lines show 
the case for n = 0.1 and 1 cm~^, respectively. For the CS dust 
echo, the solid line shows the case for M = 10~^ Mq yr~^ and 
*wind = 10^ yr (outer radius R2 = 3 X 10^^ cm). The expected 
luminosities of the SN dust and the shocked CS dust are indicated 
with the arrows. 



(LiR - 1.5 X lO^^crg s-i) of SN 1978K. 

A possible additional heating source is UV or X-ray 
emission from the cjccta-CSM interacting region. X-ray 
luminosity of SN 1978K is very high (3 x lO^^erg s"^ 
(jSmith et al.ll2007|) . being enough to account for the ob- 
served IR luminosity. However, the X-ray emission from 
the shocked region is absorbed by the optically-thick 
ejecta via photo-absorption e ven at the transitional phase 
(jFransson fc Chevaliedll987f ) , and it cannot reach the in- 
ner layer and cannot heat the SN dust there effectively. 
Therefore, SN dust at the inner layers of the ejecta is 
quite unlikely to be the origin of the IR emission from 
SN 1978K. 

Note that dust formation in a cool dense s hell has been 
propo sed for some interac ting SNe (e - g-- [Pozzo et al 
20041 iSmith et al. 2008; Mattila et ahl 120081: |Fox et al 



200a ISmith et al.ll2009i: iChugail 1200911 The cool dense 



shell may be formed at the reverse-shocked region, and 
thus, the dust can be heated by the UV/X-ray radiation. 
Although it is not clear whether 1.3 x 1O~'^M0 of the 
dust can be formed in the cool dense shell, this scenario 
cannot be completely excluded as a possible origin for 
the IR emission from SN 1978K. 

5. IR EMISSION FROM NORMAL SNE 

From the upper limits of the flux for 5 non-detected 
SNe, upper limits of the emitting dust mass can be ob- 
tained. In this section, we first summarize expected 
IR emission from normal SNe at the transitional phase. 
Then, we discuss the implication of the non-detection by 
current observations. 

5.1. Summary of IR Emission at the Transitional Phase 

For CSM-rich SN 1978K, we have discussed IR emis- 
sion expected by (1) light echo from CS and IS dust, 
(2) shocked CS dust, and (3) SN dust. In this section, 
we summarize the IR emissions from normal SNe by 
these processes, presenting the expected IR luminosity 
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and typical dust temperature. 

We have shown the hght curve of the hght echo by CS 
dust for M = 10"'' Mq yr-i in Figure H Using the 
same method, we calculate the echo luminosity with a 
normal mass-loss rate M = lO"'^ Mq yr~^, and a normal 
SN luminosity Lq ^ 6 x 10^^ erg s~^. The solid blue 
line in Figure [6] shows the IR luminosity with the wind 
duration twind = 10^ yr. The echo luminosity is about 
L = 3 X 10^6 erg s"^ at t = 30 yr. 

Figure [6] also shows the luminosity by IS dust echo as 
calculated in Section 14.2.11 The light curve of the IS 
dust echo has a long timescale, and will exceed that of 
the CS dust echo after i = 5 - 10 yr for the ISM den- 
sity of n = 0.1-1.0 cm~'^. Note th at a typical interstellar 
densit y is not well constrained. iGaustad &: van BurenI 
(|1993[ ) derived a filling factor of interstellar gas with den- 
sity > 0.5 cm~'^ to be only about 15 %. iDring et"al] 
(|1996[) derived an even smaller filling factor: 8.5 % for 
the gas density > 0.01 cm""^. When the density is lower 
than 0.1 cm~'^, the IR luminosity by the echo is roughly 
proportional to the gas density since the IS dust is opti- 
cally thin. 

For both CS and IS dust cases, the typical temperature 
of IR emission can be roughly estimated at the closest 
point from the SN, i.e., at D = ct/2 = 1.4 x 10^^{t/30 yr) 
cm, which has the largest contribution to the total lu- 
minosity. Suppose that a dust grain is irradiated by a 
SN with the luminosity LgN at a distance D. For sim- 
plicity, the spectrum of the luminous source is assumed 
to be blackbody with a temperature Tg. Then, the ab- 
sorbed luminosity is given by Labs = isN/(47r£'^)7ra^ < 
Q'^'^^Ts) >, where < Q^'^^T) > is the Planck mean of 
the absorption efficiency factor. When the dust particle 
emits radiation with a temperature Tdust, the emitted 
luminosity is L^ad = ^Tra'^aT^^^^ < Q^''"(Tdust) > ■ Un- 
der the equilibrium, i.e.. Labs = L^adi we can derive the 
temperature as follows; 



Tdust ^ 140 K 



D 



1.4 X 1019 cm 



-1/3 



L 



SN 



1 



6 X lO"'^ erg s"^ 

-1/6 



1/ 



0.1 yum 



Here we used < Q'^^'fLd nst.) > /a = O .UTTl,^,, which is 
derived for T < 100 K by[Die3 pMl . We find that this 
formula is applicable up to T ~ 250 K for our purpose. 
The expected temperature for the light echo is ^ 140 K 
at t =30 yr. 

Next we consider the shocked CS dust. Figure [5] shows 
that the typical mass of the surviving, shock-heated dust 
is about 3 x IO-^Mq at i = 30 yr for M = IQ-^ Mq yr"! 
and ■ysN = 5000 km s~^. With this mass-loss rate, the 
shocked gas density is about 6 x 10^ cm~^. For the elec- 
tron temperature about 10*" K. the expected dust tem- 
pera ture is Tdust 160 K for the grain radius of a = 0.1 
fj,m (jPwek et al.ll2008l ). The total IR luminosity can be 
derived by integrating Equation ([T]) over wavelength; 



Liu, 



il/dust 1 -^pa 



< Q^'^^CLdust) > 
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Table 5 

Upper limit of Dust Mass 



SN 


Age 


Mdust (A-/q) 


Mdust (Mq) 




(yr) 


TduBt=100 K 


Tdu,t=200 K 


Amorphous carbon 


SN 1909A 


99 


< 7 X 10-^ 


< 3 X 10-^ 


SN 1917A 


90 


< 7 X 10-1 


< 3 X 10-2 


SN 1951H 


56 


< 1 X 10-1 


< 4 X 10-3 


SN 1962M 


44 


< 4 X 10-2 


< 5 X 10-* 


SN 1968D 


39 


< 4 X 10-1 


< 1 X 10-2 


Astronomical silicate 


SN 1909A 


99 


< 3 X 10-2 


< 1 X 10-3 


SN 1917A 


90 


< 3 X 10-1 


< 1 X 10-2 


SN 1951H 


56 


< 5 X 10-2 


< 2 X 10-3 


SN 1962M 


44 


< 2 X 10-2 


< 2 X 10-" 


SN 1968D 


39 


< 2 X 10-1 


< 7 X 10-3 



3.8 g cm 



50 K 



erg s 



(7) 



where, as in Equation ([6]), we adopt < Q'^'^''(Tdust) > 
/a = 0.127T^^g^. Using this equation, the total lu- 
minosity of the shocked CS dust is estimated as ~ 
6 X 10^^ erg s""'^. For a higher mass-loss rate M = 
10^'^ Mq yr"^, the shocked CS dust mass becomes 
Mdust ~ 1 X 1O~^M0, and the temperature is Tdust 270 
K. The total IR luminosity becomes 4 x 10^^ erg s^^. 

Finally, let us consider the SN dust. For the SN dust, 
we assume Lir < 10"^^ erg s""'^, which is expected by 
the radioactive decay of ^"'Ti. If the mass of SN dust 
is Mdust = O.IM0, the temperature can be Tdust < 50 
K. Even if only Mdust = 10~^Mo of dust is formed, the 
temperature is about 100 K, which is lower than the ex- 
pected temperatures for the light echo and the shocked 
CS dust. 

In summary, IR emission of normal SNe at the transi- 
tional phase can result from three types of mechanisms, 
i.e., IS dust echo, shocked CS dust, and SN dust. IS dust 
echo and shocked CS dust have a relatively high temper- 
ature (Tdust = 150 — 270 K) while the SN dust has a lower 
temperature (Tdust < 50 — 100 K). The total luminosity of 
the IS echo can be Lir > 10^'' erg s~^ for the IS density 
of ri > 0.1 cm-^ (Figure 0. The shocked CS dust could 
also have Tir - 6 x 10^^ - 4 x 10^^ erg s-^ for M = 
lO-^-lO"" Mq yr-i (see the arrow in Figure [e]). The 
luminosity of the SN dust must be Tir < 10'^^ erg s~^. 

5.2. Constraints from Non Detections 

Figure [7] shows SEDs of 5 SNe in the 3 different host 
galaxies. These upper limits are compared with the 
model spectra of amorphous carbon (red) and astronom- 
ical sihcate (blue) with Tdust = 100 and 200 K (solid and 
dashed lines, respectively). The assumed dust mass is 
shown in the panels. In Table [SJ we show the upper limit 
of the dust mass for the case of amorphous carbon and 
astronomical silicate for the temperature Tdust =100 and 
200 K. 

The upper limit of the dust mass strongly depends on 
the dust temperature. The dashed lines in Figure[S]show 
the upper limits of the dust mass for SN 1962M as a func- 
tion of the dust temperature for the case of astronomical 
silicate (blue and green dashed lines). Constraints on the 



Observational limit for 
(astronomical silicate) 




SN ig62M 



150 
Tdust (K) 



250 



300 



Figure 8. Upper limit of the dust mass as a function of the dust 
temperature for the case of astronomical silicate. The dashed blue 
line shows the upper limits for SN 1962M, which are derived from 
the non-detection in our current observation. Solid curves show the 
upper limit for a given total luminosity L = 103^ — 103^ erg s-i. 
Our current observations are sensitive to the luminosity L > 

103*erg s-i. The observations of SNe 1 978K (pink diamond) and 
1987A (black box, IBouchet et"an[200l : IMatsuura eraI|[20TM are 
also shown for comparison. The expected ranges for the shocked CS 
dust with M = 10-^ Mq yr'^ (lower left) and M = IQ-* Mq yr-l 
(upper right) are shown (for dust size a = 0.1 — 0.3 fim). 

other 4 SNe are weaker than those for SN 1962M by a 
factor of 2-20. For comparison. Figure [S] also shows the 
estimated dust mass and temperature of SN 1978K (this 
paper) and both wari n and cold comp onent of SN 1987A 
i|Bouchet et~aI1l2004l: IMatsuura et al.ll2011[) . 

The solid lines in Figure [8] show the dust mass as a 
function of dust temperature for Tir = 10^^, 10^'', 10^^, 
and 10"^^ erg (from bottom to top). This can be 
analytically derived from Equation 

TiR 



Mdust = 5.6 X 10- 



T 



2037 Qj-g g- 



dust 

50 K 



3.8 g cm 



-3 



Mr. 



(8) 



For these lines in Figure El we do not use the analytic 
formulae, but numerically integrate Equation ([T]) over 
wavelength. 

As shown in Figure [SJ our current observations are sen- 
sitive only to the IR luminosity as high as > lO'^^erg s~^. 
This limiting luminosity is much higher than the predic- 
tions for the SN dust. The luminosity of the IS dust echo 
can be 1 X 10^^ erg s~^ with 1 cm-'^. Our current ob- 
servations can only exclude a homogeneous IS density of 
n ~ 10 cm"'^. 

For the shocked CS dust, the expected ranges of the 
dust mass and temperature are shown by two gray boxes 
in Figure [H The left lower box shows the range for 
M = 10~^ Mq yr-\ while the right top for M = 
10~^ Mq yr-^, Our current best observational limit ex- 
cludes the existence of the shocked CS dust with AI = 
10^4 Mq yr-i. 

6. FUTURE PROSPECTS 
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Figure 9. Expected IR SEDs of SN dust (purple), shocked CS dust (blue), and IS dust (red) compared with 5cr detection limit of 1 
hr imaging observation with various telescopes, i.e., JWST (black), SPICA (red), Hcrschel (black) and ALMA (black). Gray dashed line 
shows the expected confusion limit of background galaxies with a 3 m telescope (T. T. Takcuchi et al. 2011, in preparation). The distance 
to the SN is set to be 5 Mpc. For the SN dust, we assume M^ugt = O.IMq with Ttjust= 50 K. This roughly corresponds to the total IR 
luminosity of Lir ~ 10'"^ erg s~^. For the shocked CS dust, we assume 

A^dust = 3 X 1O-''M0 with 

Tdust = 150 K (the total luminosity is 

Lib. = 3.9 x 10^^ erg s~^). For the IS dust echo, we assume Lir ~ 10^^ erg s~^ with T^ust = 150 K. The effective emitting mass of dust 
is Mdust ~ 10-" M0. 



We show that our current observations are sensitive 
to the luminosity L > 10'^* erg s~^. In this section, 
we discuss prospects for future observations. Figure 
[9] shows the expected flux for the SN dust (purple, 
-^dust = O.IM0 and Tdust = 50 K), the shock-heated 
CS dust (blue, Mdust = 3 x lO""^ Mq, and Tdust = 150 
K, expected for M = 10"^ Mq yr'^), and the IS dust 
echo (red, L = 10^'' erg and Tdust ~ 150K, the effec- 
tive emitting dust mass is Mdust ^ lO^^'Afo). The SN is 
assumed to be located at 5 Mpc. 

The model spectra are compared with 5cr detection 
limit of 1 hr imaging observations with various tele- 
scopes, i.e., with The James Webb Space Telescop e 
(JWST, black, ' ' 
SPICA (red. 



Wright et all[200l lSwinv ard et al. l [200l 
Nakagawal l2'oTor Wada fc Katazal 12010: 



Swinvard '200S; Kataza et al.' '2010) He r schel (black 
Pilbratt ct al, 2010: Pogli tsch et al.i 120101: iGriffin et al.l 
2010D, and ALMA (black, iWootten fc ThompsonI 1200911 



For SPICA, we also show the limit for spe ctroscopic ob- 
serva tions (i? = 50) in the red dashed line ([Kataza et alj 
|2010[) . In Appendix |^ we show similar plots for other 
dust species. 

The figure shows that future observations at < 30/im 
are sensitive to the shocked CS dust or IS dust echo. 
These two cases may show similar overall spectra, but 
ionic emission lines are also expected for the shocked CS 
dust. Low resolution spectroscopy will be possible with 
SPICA, which will be useful to distinguish the shocked 
CS dust emission from IS dust echo. 

IR emission by SN dust has a smaller contribution to 
the IR luminosity. However, thanks to the lower tem- 
perature than the shocked CS dust and IS dust echo, 
the expected spectrum has a peak at longer wavelengths. 
The contamination by the shocked CS dust may be small 
if M < 10^^ M© yr^^. In addition, the emission from 



the IS dust echo can also be sma ll if a typical interstel- 
lar d e nsity is less than 0.1 cm"^ (jCaustad fc van Burenl 
ring et aLlll996f ). If the contributions from these 
two components are small, the detection of SN dust is 
limited only by the confusion limit of background galax- 
ies, which is an order of 0.1 mJy at 50 fim with a 3m 
telescope. 

The great advantage to detect the SN dust at the tran- 
sitional phase is the low optical depth. If a SN forms 0.1 
Mq of dust, the optical depth becomes 



— K^ Mdust ^ 

in 



iO.l 



562 cm2 g-i J V O.IMq 



dust 



VSN 



5000 km s- 



30 yr 



(9) 



For the mass absorption coefficient of dust n^, we take 
the va lue for astronomical silicate at 24 ^m (|Drainel 
l2003bD . At longer wavelengths, the coefficient is lower 
(ki/ (X A~^). At early epochs, such a large mass of dust 
is optically thick even at IR wavelengths. On the con- 
trary, at the transitional phase, the optical depth can be 
smaller than unity, which enables a reliable estimate of 
the total dust mass. 

Another advantage is the small confusion with the IS 
dust swept up by the SN ejecta. The mass of the IS 
dust (A^isdust) can be approximately written as follows 
by neglecting the destruction; 

Misdust ^ -^{vsNtfnisMmpf~^ 



;2 X 10" 



5000 km s- 



30 yr^ 
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/^^isM \ /Xy^Me. (10) 
Vl cm-^J \200 J ° ^ ' 

The mass of swept- up IS dust is nearly O.IMq in SNRs 
ai t = 1000 yr, which makes it difheult to estimate the 
SN dust mass in SNRs. On the contrary, the swept-up 
IS dust mass at the transitional phase is only an order of 
IG-^Mq. 

In summary, future observations of SNe at the transi- 
tional phase will the gap of infrared observations of SNe 
with the age of 10 — 100 years. At < 30 /im, observations 
will provide an interesting way to study the circumstel- 
lar and interstellar environments of the progenitor. With 
SPICA, these components can also be studied by low 
resolution spectroscopy. If these two components are rel- 
atively weak, the SN dust might be detected at longer 
wavelengths, which enables the reliable measurement of 
the dust mass. 

7. CONCLUSIONS 

SN explosion have been studied mostly by observations 
of extragalactic SNe at the early phase {t < a. few years) 
or observations of SNRs at the late phase (t > 100 yr). 
Observations at the transitional phase from SN to SNR 
have not been extensively performed, especially at IR 
wavelengths. We present theoretical predictions for the 
IR emission from SNe at the transitional phase. We show 
that the emission arises from SN dust, light echo by CS 
and IS dust, and shocked CS dust. 

We search for IR emission toward 6 core-collapse SNe 
at the transitional phase in the nearby galaxies NGC 
1313, NGC 6946, and MlOl by using the data taken with 
the AKARI satellite and Spitzer. Among the targets, we 
detect emission associated with SN 1978K in NGC 1313. 
The emission is explained by 1.3 x 10"^ A/0 of silicate 
dust. The shocked circumstellar dust for gas mass-loss 
rate M = 10~^ Mq yr~^ is the most probable origin 
of the IR emission from SN 1978K. Dust formed in the 
outer, cool dense shell might also be a possible origin. 

IR emission from the other 5 objects is not detected. 
Current observations are sensitive only to the total lu- 
minosity of L > 10^® erg s~^. The non-detection of SN 
1962M excludes the existence of the shocked CS dust for 
the high gas mass-loss rate AI = 10~^ Mq yr~^. 

Future observations will fill the gap of the IR observa- 
tions at the transitional phase. At < 30 fim, the emission 
is dominated by the shocked CS dust and IS dust echo. 
The observations will provide an interesting opportunity 
to study the circumstellar and interstellar environments 
of the SN progenitor. SN dust component has a smaller 
contribution, but has a spectrum peaking at longer wave- 
lengths. If the emission by the shocked CS dust and IS 
dust echo does not dominate the IR emission, the SN 
dust might be detected at > 30 /xm. This will give a reli- 
able estimate of the SN dust mass, which does not suffer 
from the high optical depth of dust and the confusion 
with the interstellar dust. 
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APPENDIX 

SEDS FOR AMORPHOUS CARBON, GRAPHITE, AND FORSTERITE 

As discussed in Section |6l future mid/far-IR observations will be able to detect the shocked CS dust, IS dust echo, 
and possibly SN dust. Figure |9] only shows the case for astronomical silicate. In Figure [TOl we show similar plots, but 
for forsterite, amorphous carbon, and graphite. For each component, the same temperature and dust mass with those 
for astronomical silicate are assumed. 
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Figure 10. Expected SEDs of SNe at the transitional phase [t = 30 yr) for forsterite, amorphous carbon, and graphite (from top to 
bottom). For the shocked CS dust, Tju^t = 100 - 150 K and M^^^t = 3 X IQ-^Mq are assumed. For the IS dust echo, T^j^g^ = 150 K and 
-'^dust = 1 X 1O~*M0 are assumed, which give the luminosity of L = 10^''' erg s~^. For the SN dust, T^^at = 50 K and M^uat = O.IMq 
arc assumed. 



